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Abstract 

We present a general theory for the equilibrium current distribution in an 
interacting two-dimensional electron gas subjected to a perpendicular mag- 
netic field, and confined by a potential that varies slowly on the scale of the 
magnetic length. The distribution is found to consist of strips or channels 
of current, which alternate in direction, and which have universal integrated 
strength. 
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1. Introduction 

The quantum Hall effect [|I|] is described as the quantization of the off-diagonal conduc- 
tance of a two-dimensional electronic device, at integral or simple fractional multiples of 
e 2 /h, and the vanishing of the diagonal conductance, when a perpendicular magnetic field 
of appropriate strength is present. It has long been understood |2j that this "appropriate 
strength" of the magnetic field is such that the the Fermi level lies in a gap of the spectrum of 
bulk extended states. Finer features of the effect are not nearly as well understood. There is 
considerable confusion, for instance, concerning the question of the physical distribution of 
the Hall current in a quantum Hall device. Is the Hall current carried by "edge states" resid- 
ing near the the edges of the device, or is it carried by the bulk of the electron fluid? While 
early work ||, based on a model of non- interacting electrons, supported the first alternative, 
more recent work [|j] taking into account Coulomb interactions self-consistently, concludes 
that at least a substantial part of the current is carried in the bulk. Experimental attempts 
to resolve the ambiguity by directly imaging the voltage and/or current distribution, have 
been made with no conclusive results so far. A major practical difficulty is that in a 
nonuniform electron fluid subjected to a magnetic field, a finite current density exists even 
in thermal equilibrium, due to the time-reversal symmetry breaking action of the magnetic 
field. Although the net current through the device is zero, its local value may be quite large 
- in fact larger than the Hall current density, defined as the difference between the total 
current density measured in the transport regime and the current density at equilibrium. It 
is evident that a clear understanding of the properties of the equilibrium current distribution 
is a prerequisite for the understanding of experiments on the Hall current distribution. In 
two recent papers ||[7j we studied in detail the properties of this distribution, and found 
that they could be summarized in a few fairly general statements. We considered a confining 
potential that varies slowly on the scale of the magnetic length I: V'l « A where A is 
the minimum energy gap that we want to to take into account (A = hu c if the fractional 
quantum Hall effect is ignored). This condition - at least insofar as the integral Hall effect 
is concerneed - is very well satisfied in realistic devices, in which the confinement lengths 
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are of the order of 10 4 A or more, while I is of the order of 100 A |§. Our expression, which 
correctly describes the current distribution on a scale larger than /, has two components: an 
"edge" current, which is proportional to the density gradient, and a "bulk" current, which is 
proportional to the gradient of a self-consistent potential. In the limit of zero temperature 
these two contributions become separate in space: The edge current flows in the "com- 
pressible" regions, in which the density varies, while the self-consistent potential remains 
constant. The bulk current flows in the "incompressible" regions, where the self-consistent 
potential varies and the density is constant. The directions of these two types of currents 
generally display a striking alternating pattern. Furthermore, the integrated currents across 
a compressible or incompressible region are universal, in the sense of depending only on the 
chemical potentials of a uniform two-dimensional electron gas near an incompressible state. 

2. Current Density Distribution 

A convenient tool for the study of equilibrium orbital currents in an interacting electronic 
system is provided by the current-density functional theory (CDFT) ||. This formulation 
provides a rigorous mapping of the many-body problem for the ground-state energy, density, 
and current density, to an effective one-body problem, in which independent electrons move 
under the action of self-consistent local scalar and vector potentials. Unlike the ordinary 
density functional theory, CDFT is rigorous, in principle, for the calculation of equilibrium 
currents. The effective single particle equations have the form || 

( -L[-iV + -A{r)f + V{r) + V H {r) + V xc {r) + [-zV + -A[f)\ ■ ^A xc {r) 
2m c c 2mc 

+ ^-A xc {r) ■ [-zV + -A(r)] - f -K0 ' ^ = 
2mc c c n{r) 

where A(r) and V(f) are the external vector and scalar potentials, n(r) is the number 
density, and j (f) is the current density. Here 




is the Hartree potential, 
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_ 8E xc [n{^B v {r)] 
Vxc[r} ~ 6n(fj [6) 

is the exchange-correlation (xc) scalar potential, and 

e_ A = _mc f 6E xc [n(f),B„(r)] 
c V ; en(r) 8B v {f) 

is the exchange-correlation vector potential. A crucial feature of the theory is that the xc 

energy functional E xc [n(f), B v (r)] is a functional of the density n(r) and the vorticity || 

4(0 = S OT -™^xfM), (5) 
e \n(r) / 

where B(r) = V x A(f) is the magnetic field. The density and orbital current density 
are self-consistently determined by sums of orbitals, which are solutions of eq. We 
have ignored spin for simplicity. The theory would yield the exact density and current 
distributions, if the exact exchange-correlation energy functional were known. In the local 
density approximation (LDA), which is justified for density distributions that are slowly 
varying on the scale of the magnetic length, E xc takes the form 

E xc [n{f),B u {f)} ~ J e xc [n(r),B u (r)}df (6) 

where e xc [n, B] is the xc energy density of the two-dimensional electron gas (2DEG) of 
uniform density n in a uniform magnetic field B. 

We now specialize to the case of a uniform magnetic field perpendicular to the plane of 
the electrons A(f) = Bxy in the Landau gauge, and a slowly varying confining potential 
V(f). Following the approach of |J, at each point r in space we define local cartesian axes 
such that the x axis is along the direction of the gradient of the total self-consistent potential 
Vhxc = V + Vh + V xc . Thus, the total potential is locally a function of the x coordinate only. 
Choosing also the Landau gauge in such a way that A depends only on the local x coordinate 
(and is parallel to the local y axis) we effectively obtain a local one-dimensional problem. 
In particular, the density and current profiles depend only on the local x coordinate. This 
is so, because an x-dependent density and current profile leads to an x-dependent xc vector 
potential in the y direction given by the formula 



e _ A ( x \ _ mc d SE xc [n(x),B u (x)} _ 
c xc en(x) dx 5B p (x) 

It is then clear that the eigenfunctions of the Kohn-Sham equation localized near r*o have 
the form 

^nk(x,y) = ^ r e iky ij nk (x), (8) 

VLy 

where the ip n k( x ) are normalized solutions of the equation 

(~l~Ix* + V X ~ X)2 + Vhxc ^ + ~ C A ^ X )\- X ~ X ~ f^l) ^k{x) = e nX ^ nk (x), (9) 

where x is in units of /, n is in units of l~ 2 , j y is in units of h/2ml 3 , Vh xc is in units of 
hu c , eA xc /c is in units of h/l, and X = —kl is the guiding center of the Landau orbital). A 
more rigorous discussion can be given in terms of exponentially localized Green's functions 
0. Eq. (^|) may be solved perturbatively, which is exact in the limit of slowly varying 
potential. Because the eigenfunctions of —d 2 /dx 2 + (x — X) 2 are strongly localized about 
x = X we can expand 

V Hxc {x) « V Hxc (X) + {x- X)Vti xc (X) (10) 

and use 

-A xc (x)[x -X- « -A xc (X).[x - X] (11) 

c n[x) c 

In writing down the second equation we have exploited two facts: (i) A xc (x) being given by 
eq.(|7|) is already a gradient of a local function of x. Hence, we truncate its expansion to the 
zeroth order in x — X (ii) j y {x), as it will soon be clear, is also a gradient of a local function 
of x, therefore its multiplication by A xc (x) is the product of two gradients, which must be 
regarded as a higher order infinitesimal quantity and hence dropped in the present theory. 

Using equations ( [T0| ) and (0), equation (||) can be solved, and the density and current 
density can be computed explicitely. The general character of the density distribution has 
been discussed at length elsewhere P|||7|]. Very briefly, the self-consistent density profile 



consists of alternating compressible and incompressible regions. A typical profile is shown 
by the solid line in Fig. 1. In a compressible region, the density varies while the self- 
consistent potential is pinned to a constant value. In an incompressible region, the density 
is pinned to an integral or fractional multiple of l/27r/ 2 , and the self-consistent potential 
varies. The self-consistent chemical potential is tied to a set of quasi-degenerate single 
particle orbitals, which are localized in the compressible regions. Rather than insisting on 
the properties of the density distribution, which are well understood, here we focus on the 
result for the current density, which can be stated as follows: 

dTl 6 

3v(z) = 7fc^ + 2n(x)[V^ xc (x) + -AJx)\ (12) 

where 

7fc [n,S] = (2M + l), (13) 

and \y\ is the integral part of the filling factor. As a final step, we observe that from eq. 
(0), combined with the local density approximation, eq. (0), we obtain 

2-A xc [x)n[x) = 7TC [n,B]^, (14) 

where we have defined 

1 <9/i xc [n, B] 

7«M = -— ^— . (15) 

Here fis = eh/2mc is the Bohr magneton and \x xc = de xc [n, B]/dn is the xc part of the 
chemical potential of the uniform 2DEG. We have set B v = B after taking the functional 
derivative of E xc , since A xc is already of the desired order in the gradient of the density. 
It is important that, in the definition of ~f xc [n,B], the density is never taken to be exactly 
equal to one of the incompressible densities at which \h xc is discontinuous. This prescription 
assures that ^ xc is finite everywhere, although possibly discontinuous at the incompressible 
densities. 

Restoring the physical units and returning to the original coordinate system we finally 
write our result as 
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3(f) = 7?-{lk[n, B] + >y xc [n, B])Vn(r) x z + -^—n(r)VV Hx c(r) x z. (16) 

To understand this formula, we must keep in mind that, in the low temperature limit, the 
"bulk" term (W(r)) and the "edge" term (Vn(r)) contribute to mutually exclusive regions 
of space - incompressible and compressible regions, respectively. In the incompressible 
regions, the density is pinned to one of the quantized values, and the bulk term reduces 
to the usual Hall current j = (ei>i/h)WHxc(f) X z, where z/j is a quantized filling factor, 
and Whxc plays the role of the electric field. Assuming, for the sake of definiteness, that 
the self-consistent potential increases monotonically from the center to the perifery of a 
quantum Hall system (for example a bar or a dot), this contribution to the current is found 
to be paramagnetic. In the compressible regions, the self-consistent potential is pinned, 
while the density varies. Here the edge term becomes operative. Under the assumption 
that the density decreases monotonically from the center to the perifery of the system, its 
contribution to the current is found to be diamagnetic. A typical self-consistent current 
density profile is shown ClS db dashed line in Fig. 1. 
3. Orbital Magnetization 

Since the charge current J(f) = —ej(r) satisfies the continuity equation 

V ■ J(r) = 0, (17) 

it is possible to write it as the curl of a local function M(f)z which we call the "orbital 
magnetization" : 

J(f) = cV x M(r)z. (18) 

It is evident that M(f) is defined only up to an arbitrary additive constant. However, in the 
limit of slowly varying density, the orbital magnetization may be derived from the internal 
energy of the uniform electron gas 

This gives [[UJ 



with e = €k + e xc the total internal energy (kinetic plus exchange-correlation) per unit area 
of the uniform 2DEG. It is not immediately evident that eq.(|2"U|) is equivalent to eq.(|T^). A 
puzzling feature is that the VV(r) term does not appear explicitly in eq . (|20|) . The puzzle is 
resolved by noting that the thermodynamic derivative appearing in eq. (|20|) has singularities 
when the density is such that the compressibility vanishes. These singular contributions are 
responsible for the appearance of the VV(r) term in the incompressible regions. To show 
this, we rewrite d 2 e/dndB as 



d 2 e[n,B] [dpi] (dn/dB) 



(21) 



dndB \dB J (dnjd^B 
When the density is not such that the compressibility vanishes, the derivatives are finite, 
and it is easy to verify that eq. ( f20l) reduces to the first term of eq. ([il|). But when the 
density is such that the compressibility vanishes, we must turn to the second equality in 
(pl|), for insight into the behavior of the current. The denominator (<9n/<9/i)# vanishes, 
because the system is incompressible. At the same time Vn(r) also vanishes in the regions 
where incompressibility holds. The ratio of these two quantities, however, remains finite in 
the limit T — > 0: 

-(|^) Vn(r1 - VV Hxc {r). (22) 

The numerator (dn/dB) ^ is well known to be proportional to the Hall conductivity e 2 Vi/h 
|nj. This establishes the equivalence of eqs. (P0|) and fllTf). 

The usefulness of writing the current as the curl of the orbital magnetization lies in the 
fact that it makes possible to express the flux of current across an arbitrary path as the 
difference in the values of the orbital magnetization at the end points: 

h,2 = I" J(r) • n{f)dl = c[M(l) - M(2)], (23) 



where the line integral is along a path joining points 1 and 2, and n is the unit vector normal 
to the path. The key observation is that the orbital magnetization at the end points may be 

8 



well represented by the LDA even if it is not well represented at intermediate points along 
the path. 

As a simple example, consider the disk of non-interacting electrons studied recently 
by Avishai and Kohmoto [12j (AK). The integrated current across a radius is simply the 
magnetization at the center of the dot, since the magnetization obviously vanishes at the edge 
of the system. But the density is uniform at the center of the dot, so the LDA is certainly 
valid there. Furthermore, under the conditions envisaged by AK, the Fermi energy is pinned 
to the bottom of the i-th Landau level in the bulk, which means that \x = (i + l/2)hu^ 
(a superscript "+" or "-" means that an infinitesimal positive quantity must be added 
or subtracted). Looking at the relation between n and /i in the non- interacting 2DEG at 
an infinitesimal temperature we see that this chemical potential corresponds to a density 
n — [(i + l)/27r/ 2 ] - at the center of the dot. This, in turn, corresponds to an orbital 
magnetization M = — (2/zb/27t/ 2 ), as one can see from Fig. 2. Therefore the total current 
integrated from the center to the edge of the dot is equal to — ieu c /4:ii, i.e., the total current 
is quantized in integral multiples of eu c /47r. This result was first derived by AK from a 
numerical solution of the Schrodinger equation with hard wall boundary conditions. The 
present derivation demonstrates the generality of the result. 

Next we consider the interacting electron fluid. The energy of the uniform phase consists 
of kinetic and exchange-correlation contributions. The kinetic energy is known exactly, and 
the kinetic contribution to the orbital magnetization is given by 

M k (u, B) = (H - (2H + - H) ^L. (24) 

This curve is plotted in Fig. 2. No simple expression is known which represents e xc ac- 
curately at all values of filling factor v. However, to calculate the integrated current in a 
compressible or incompressible region, we only need to know what happens in the vicinity of 
an incompressible filling factor uq, integral or fractional. Here some essentially exact results 
can be obtained. In a neighborhood of uq the exchange-correlation energy density can be 
expanded as follows: 
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txcO + 



2tt/ 2 



li xc [v£,B] + 0((u - u 



,3/2N 



(25) 



for v > z/ , and 



(26) 



for v < uq. Here e xc o is the xc energy density of the incompressible state at filling factor 
i> and HxcIvq, B] and ^i xc [i/q,B] are the xc chemical potentials calculated as right or left 
derivatives of the xc energy density with respect to density. In the limit of high magnetic 
field (e 2 /lhu c « 1) e xc0 scales as e xc0 (B) = e xc0 e 2 /27rl 3 and the /Ws scale as fi xc (u,B) = 
fixc(y)e 2 /I, where the dimensionless quantities e xc0 and Jjl xc {u) are independent of magnetic 
field. The 0((u — z/o) 3//2 ) term on the right hand side represents the interaction energy of a 
classical Wigner crystal of dilute quasiparticles (or quasiholes) added to the incompressible 
state. It is easy to verify that this term - as well as all the higher order terms in the 
expansion of the energy in v — u - does not contribute to the orbital magnetization in the 
limit v — > u . Taking the derivatives of eq. (p5|) and eq. ( p6|) with respect to B, using the 
power law scaling of the energy and chemical potentials as functions of B, we obtain the 
exchange- correlation contribution to the magnetization 



M xc (v = p ,B)= -3e xc0 + 1vqii xc {v q 



(27) 



lhuj c 27r/ 2 

The above equation may be used to calculate the exact total current that flows through 
a compressible strip or channel connecting an incompressible region at filling factor u to 
another incompressible region at filling factor V\. Assuming, for definiteness, that uq > v\ 
and integrating the current flow along a line going from the higher to the lower density we 
find that the total current is 



I = c \M{v = Uq, B) - M{y = u+, B)), 



(28) 



where M = + M xc . The remarkable feature of this result is that it holds true despite the 
lack of detailed information about the density and current distribution within the compress- 
ible region, i.e., despite the fact that n{f) and 7 xc (n, B) are not known everywhere inside 
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the region. Equation fl27j ) can also be used to calculate the current flowing through the 
incompressible region, say at filling factor V\ . For this region we obtain 



where Ai = fi(h>^~,B) — fi(ui,B) is the gap in the electron chemical potential at filling 
factor v\. This result is precisely what one would expect, based on the observation that 
the change in the self-consistent potential across an incompressible strip must equal the 
incompressibility gap. 
4. Conclusions 

In conclusion, we have elucidated in this paper the nature of the equilibrium current dis- 
tribution in a two-dimensional electron gas confined by a potential which varies slowly on the 
scale of the magnetic length. We have shown that the density profile consists of alternating 
compressible and incompressible regions, and that the current flows in opposite directions 
in these two types of regions, reflecting the sawtooth behavior of the magnetization of the 
uniform interacting electron gas as a function of density. Finally, we have demonstrated that 
the integrated current across a compressible or incompressible region can be expressed (in 
the limit of strong magnetic field) in terms of the chemical potentials of the uniform electron 
gas near incompressible filling factors. Detailed numerical results for exchange-correlation 
chemical potentials near integral filling factors will be reported elesewhere. 
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FIGURES 

FIG. 1. Typical self-consistent equilibrium density (solid curve) and current density (dashed 
curve) for a two dimensional electron fluid confined along the x direction at low temperature. The 
density is plotted in units of po = \/2irl 2 and the current is plotted in units of jo = eu> c /2irl. Note 
the alternating signs of the edge and bulk currents 

FIG. 2. Orbital magnetization density versus filling factor for a non- interacting uniform 
two-dimensional electron gas at low temperature. 
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